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Abstract 

The van Hove autocorrelation function of the impurity is obtained which is exact to the first 
order in the vacancy concentration. It is found that in the case of strong vacancy-impurity binding 
a singularity in the van Hove function corresponding to a resonant bound state develops on the 
unphysical sheet in the complex frequency plane close to the real axis. It is argued that this bound 
state corresponds to the defect-impurity pairs widely used in models of diffusion in semiconductors. 

PACS numbers: 66.30.Dn, 66.30.Jt 
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FIG. 1: Frequencies of exchange of the vacancy (grey square) with atoms in in the vicinity of the 
impurity (black circle) as defined in the five-frequency model: wi-exchange with the atoms in the 
first coordination sphere; ^-exchange with the impurity atom; u^-dissociative jumps away from 
the impurity into higher coordination spheres; ^-associative jumps from higher spheres into the 
first one. Not shown in the figure is the frequency ujq of exchange with the atoms in the host bulk. 



In close packed solids (like FCC crystals) the diffusion of atoms goes predominantly via 
their exchange with vacancies. Because the vacancy concentration is usually small, in the 
limit of infinite dilution a microscopic description of the diffusion reduces to the consideration 
of the two-body interactions and exchanges between a single vacancy and the atom of the 
diffusing specie 

When microscopic interactions are known, a two-body problem can usually be solved 
exactly. For example, in Ref. a solution for the tracer Green function accurate to the 
leading orders in impurity and vacancy concentrations was obtained in the case of the self- 
diffusion problem. This success was due to the simplicity of the self-diffusion case where all 
vacancy hopping parameters are identical. 

In the case of impurity, however, the vacancy in the host bulk and in its vicinity expe- 
riences different forces, so the number of parameters grows. This case is most frequently 
described in the framework of the so-called five-frequency model (5FM; see Fig. |TJ) []], 0, [| . 
Introduced by Lidiard [^| half a century ago, this model has became a conventional tool for 
the microscopic description of impurity diffusion. 

The aim of the present Letter is to obtain the solution of the 5FM of impurity diffusion 

n 

similar to that of Ref. |4| for the tracer diffusion and to show that under strong vacancy- 
impurity(V-I) binding the solution contains a resonant bound state which can be identified 
with the defect-impurity pairs widely used in the description of diffusion in semiconductors 



BQ. 
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To avoid awkward matrix formalism conventionally used in microscopic theories of diffu- 
sion , the calculations of the present Letter will be based on the second quantization 
representation of stochastic dynamics 8]. This allows for the treatment of diffusion prob- 
lems within the standard techniques of the quantum theory. A detailed exposition of this 
approach may be found in the review article jjj. The formalism which will be used in 
this Letter is presented in detail in Ref. so below are listed only the formulas used in 
subsequent calculations. 

In the second quantization approach the host lattice completely filled with the solvent 
atoms is considered to be the ground state (|host)) of the diffusion problem. The impurity 
atoms and the vacancies are created (annihilated) on lattice site % under the action of the 
so-called hard-core boson creation (annihilation) operators af (<0 and bf (6j), respectively. 
These operators satisfy the canonical bosonic commutation relations in all cases except when 



on the same site in which case they satisfy the anticommutation relations 



11] 



{4,di} = l and d 2 = (d+) 2 = (1) 

(here di stands either for a, or for bj) to prevent multiple occupancy of the site by the same 
specie. The particle number operators as usual are 



n 



afdi and Vi = bfbi (2) 



and the average concentrations are their statistical averages Cj = (n*) and cy = (vi). 
The coherent state 

i)=n(i+^)(i+^)i h ^t), (3) 

i 

comprises all possible configurations of impurities and vacancies with equal weight and is 
used below to compute statistical averages. The weak equivalences between the operators 
with respect to the conjugate coherent state 

<k = (k; flat = <|(i-70 (4) 

and similar relations for the vacancy operators are easily derived from Eqs. (jT))-(J3{). 

The matrix of transition probabilities T governing the time evolution of a state vector 

|t)=exp(Tt)|0>, (5) 
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obeys the condition of probability conservation 
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(|T = 0. (6) 

The van Hove autocorrelation function Giq describes the evolution at equilibrium of an 
impurity atom initially placed at the origin of the lattice coordinates. It can be obtained as 
the Cj — > limit of the impurity- impurity correlation function because in this limit only the 
autocorrelations survive: 



Giq = lim 

cj— »0 



J - 1 <hexp(tTHe- w / fer |)l, (7) 



12j, the transition matrix 



where 7i is the Hamiltonian of interparticle interactions which includes a normalization 
constant such that (\e~ n / kT \) = 1 [10]. With the use of Eq. © the last equation can be re- 
written as a conventional retarded correlation function. For brevity, the following notation 
for such functions will be used: 

(A; B) t = hm [cj 1 (\A(t)B(0)e~^ kT \)] , (8) 

where A(t) = e~ tT Ae tT . Thus, the evolution of the correlation function can be studied with 
the use of the Heisenberg equation 

A = [A, T] (9) 

provided the transition matrix is known. According to Ref. 
consists of two parts: 

T = T in -T out , (10) 

where the first part describes intersite particle exchange while the second part is site-diagonal 
and depends only on the operators rij and v j. In the case under consideration the "in" part 
can be written as 

I* = EK+i^ + a,)^. (11) 

Here the first term in the parenthesis corresponds to the vacancy exchanges with the host 
atoms while the second term describes the V-I exchanges. In general case the hopping 
frequencies entering Eq. (jllj) depend on the whole many-body configuration of vacancies 
and impurities. But in the double dilute limit under consideration it is sufficient to keep 
the dependence of u>ij on the position of only one impurity, while ufj does not depend on 
the configuration at all because the two participating particles are already present in the 



transition matrix. So in the two-body approximation ufj can be identified with 002 while uiij 
have the values shown in Fig. It is important to note that Uij 7^ ujji. 

The "out" term in Eq. (jlUj) can be easily obtained from Eq. (JHJ) with the use of Eq. 
(jlj). But because it commutes with all site-diagonal operators, it will not be needed in the 
calculations below. Thus, substituting A = and T m into Eq. (JSJ) and computing the 
average in Eq. (JHJ) with B = n to the leading orders in ci and cy, with the use of Eq. (HJ) 
it is straightforward to obtain 

dG m (t)/dt = w 2 E (( n i+sW, n ) t - (nivi+s] n )t) , (12) 
s 

where by 8 are denoted 12 different lattice vectors connecting nearest neighbor sites. 

Similarly, by choosing A = riiVj and B = n Q in Eqs. (jHI) and Q one arrives at a closed 
set of equations necessary to calculate the impurity van Hove function Gio(t) to the chosen 
order of approximation 

(d/dt){niVj; n ) t = ^(rijVi; n ) t - ufjiniVf, n ) t 

+ ^2^ji(niVi;n ) t - (^^(niVf, n ) t . (13) 

i i 

The above two equations have the same formal structure as the corresponding equations of 
Ref. 0|, so their solution can be found by a straightforward use of the techniques of that 
reference. Therefore, only notation and the points of departure from the solution of Ref. 4] 
will be explained in some detail. 

The first step is to subject the equations to the Laplace (with respect to t) and spatial 
Fourier transforms (L-F transform). Their normalizations should be understood from the 
definition of f K (r) (cf. Ref. |j) 

(n,tw; no), = ^ / ^ E ( 14 ) 

and a similar expression for Gio(t)-Gic(z) pair. The Laplace transform of derivatives contains 
the function values at t — 0. In the present case they are given by the equilibrium correlators 
in (JHJ which are G i(0) = 5 10 and 

(mvjnoh = c v 5 m {l - <%) [1 + (e Eb/kT - l),5\i- jm ], (15) 

where the multiplier in the parentheses forbids simultaneous occupation of the same site by 
the impurity and the vacancy while the term in square brackets describes the enhancement 



(suppression) of the vacancy concentration in the first coordination sphere of the impurity 
provided the V-I binding energy Ef, is positive (negative). The transformed Eq. (|T2*j) reads 
(cf. Eq. (5) of Ref. Q) 

zG K (z) = 1 - oo 2 Y\l - exp(tK ■ 5)]f K (d). (16) 

According to Ref. it should be cast into the form of the diffusion propagator 

G K (z) = [z + c v F(z,K)e K ]-\ (17) 

with the z and l^-dependent diffusion correlation factor 

F(z, K) = ^Y\ l ~ *MiK ■ S)]f K (S), (18) 
Cv^K s 

where €k = ^oJ2s [1 ~~ exp(iK ■ 5)]. 

Eq. (llHjl under the L-F transform takes the form 

zf K {r) = c v [l - S(r)][l + (e& _ i) £ S (r - d)} 

s 

+u 2 J2$(r - S)[exp(-iK ■ 6)f K {-r) - f K (r)] 
s 

+ ^ Mr - 6, 6)f K (r - d) - w (r, S)f K (r)} , (19) 
s 

where u(r, 8) is the frequency of hopping from site r to r + d. It is easy to see that 
/k(0) = satisfies this equation irrespective of /kt(t*) values at other sites because in 
the 5FM the vacancy hops onto or from the site occupied by the impurity are forbidden: 
u(—8,6) = u(0,6) = 0. Therefore, henceforth /j<-(0) will be excluded from the considera- 
tion. Furthermore, because in the absence of fluxes (K = 0) the system is at equilibrium, 
the values of fK=o( r ) should be equal to the L-F transformed Eq. (fT5|). The result can be 
written symbolically as 

fo = E/z, (20) 

where the vector components are the values of corresponding quantities at lattice sites r. 
Because the components of E coincide with the first line on the right hand side (RHS) of 
Eq. (fl9*|) . zf turns this line into identity. Thus^in order to satisfy the rest of equation the 
frequencies should be subject to the constraint |l|, |2j, |^ 

iVz/u^ = exp(—E b /kT) = e (21) 



as can be checked by a direct substitution of /o from Eq. (j2U|) into the last line of Eq. (JT§j) . 

The method of solution for general K proposed in Ref. |4J essentially consists in subtrac- 
tion from both sides of Eq. of the term which formally coincides with the equation last 
line but with all to(r, r) replaced with loq. The result can be symbolically written as 

R- 1 f K = E + Vf K , (22) 

where matrix R is composed of matrix elements 

1 expf— i\ • (r — r )] , , 

R ry = N$ ;+L ■ (23) 

the first line of Eq. ()19j) is represented by E and the remaining terms are gathered into 
V/k- The meaning of this transformation is that now on the RHS of Eq. ()22|) only a finite 
number of fK{ r ) (54 in this case) with r belonging to the first four coordination spheres 
remain because beyond these spheres u>(r, 6) — ujq = 0, as can be seen from the definition of 
frequencies in Fig. [T] Now multiplying both sides of Eq. (J22j) with R 

f K = RE + RVf K (24) 

and retaining only those equations which have at their left hand sides the same /k-(t") that 
are present on the RHS one obtains a linear system of 54 equations which can be solved, e. 
g. with the use of the Cramer's rule as 

f K (r) = A r (z,K)/A(z,K). (25) 

Here A(z, K) is the determinant of the system and A r (z, K) is the same determinant with 
the column corresponding to /k(^) replaced by RE. Substituting /k(<5) thus obtained into 
Eqs. (fT8|) and ((Zj) one obtains Gk(z) for general values of z and K. It can be used, e. 
g., to study the diffusion of Mossbauer impurities where the knowledge of the van Hove 
autocorrelation function at finite values of z and K is essential jl^l ]. 

If one is interested only in the diffusion limit z, K — > 0, then the size of the system 
(|2*4^1 can be reduced to 13 by the choice of a high symmetry direction K = (K, 0, 0) 0. 
It can be shown that the correlation factor E(0, 0) in Eq. (jl8)l can be expressed through 
the ratio of two determinants of size 13. The corresponding expression was derived and 
numerically checked on thousands of randomly generated frequency quintets. In all cases 
excellent agreement with the approximate expression due to Manning [2| was found. In 



particular, it was numerically confirmed that in the case of strong V-I binding the correlation 
factor is always enhanced as 0(l/us). This might have been guessed on the basis of Eq. (J2UJ1 
which states that /k=o(^) — £ ~ l I ' z ^ 1/^3- A careful analysis shows that in Eq. (J25)) the tu^ 1 
factor can originate only from the denominator because in the numerator the contribution 
due to RE containing e~ l is suppressed in the limit z — > by the factor z in Eq. (|18|) (see 
Eq. f)27j) below). Thus, in the case of strong V-I binding the determinant A(0, 0) should be 
proportional to u 3 . Then it is reasonable to assume that at small z and \K\ 

A(z, K) « de + C 2 z + C 3 K 2 , (26) 

where Cj. are some constants. From here it would follow that F(z, K) develop a pole-like 
singularity at some small value of z. 

To qualitatively assess the influence of this singularity on the impurity diffusion let us 
consider the simplest model with V-I binding: the so-called 2-frequency model (2FM) 

BBS 

where all frequencies are equal to ujq except uj 3 which according to Eq. (}2"T]) should be equal 
to eujq. In this case the dimension of determinants reduces to 3 so with the use of exact 
relations between R r r > derived in Ref. w| one arrives after some algebra to an expression 
which to leading orders in the small quantities reads 

F 2FM (z, K) « y 27 
z + 2ge+ (u /6)(aKy 

where g = (2AR ^)~ 1 - Because R 0: s has a square root singularity at z = jl4||, the pole 
shifts from the real axis on unphysical sheet on the distance 0(e 3 ^ 2 ) <C 1, as can be seen 
from the denominator of Eq. (|27|) . 

The quasiparticle content of the van Hove function can be established by the substitution 
of F 2 fm into Eq. (fT%|) and by approximating g with a real constant go — g{z — 0). In this 
case the van Hove function acquires simple two-pole structure 

Gk(z) „ i^w + (28) 

where 

z 1 (K)^^u (aK) 2 = D I K 2 , (29) 

z 2 (K) ^2g E + (LU /Q)(aK) 2 , (30) 

and 

4>{K) = (c v /12e)uj 2 (aK) 4 /z 2 2 (K). (31) 




FIG. 2: (color online) One-dimensional diffusion along (1,0,0) direction of the FCC lattice of initial 
5- function impurity distribution calculated in the two frequency model (see the text) for different 
values of diffusion length = 2(Djt) 1 ^ 2 . For definiteness, the parameters corresponding to the 
arsenic impurity in silicon at temperature T = 800 °C were used: Ej, ~ 1.2 eV jla ]. cy = 10 -8 (if] ], 
and a = 5.43 A. 

Thus, there are two diffusion poles. The one at —z\ corresponds to the interaction-enhanced 
impurity diffusion and the pole at —Z2 can be identified with the V-I pair state. Indeed, being 
the part of the impurity autocorrelation function this state obviously contains an impurity. 
On the other hand, its diffusivity is independent of c v which means that the vacancy is 
always present in this state. This is further confirmed by the diffusion profiles computed 
as inverse L-F transform of Eq. ()28j) and shown in Fig. 
exhibit characteristic tails which in pair diffusion models 
defect-impurity pairs. 

Unfortunately, these profiles are not quite physical because at large distances they acquire 
small negative values. More accurate calculation of the inverse L-F transform is needed to 
deside on whether this deficiency is because of too crude approximations made above or that 
0(c v ) approximation is insufficient at large distances. 
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At small diffusion lengths they 

Q 

, 121] are attributed to the diffusing 
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the work. 
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